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We present the results of searches for B decays to charmless two-body final states containing r/' 

or LO mesons, based on 20.7 fb~^ of data collected with the BABAR detector. Wc find the branching 
fractions B{B+ ri'K+) = (70 ± 8 ± 5) x 10-^ B{B'^ v'K^) = (42ti? ± 4) x 10-^ and 
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B{B+ LJTT+) = (e.eli.g ± 0.7) X 10"'' where the first error quoted is statistical and the second 
systematic. We give measurements of four additional modes for which the 90% confidence level upper 
Hmits are B{B" luK°) < 13 x 10'^, B{B+ 77^+) < 12 x 10"^ iuK+) < 4 x 10-^ 

and B{B° lotv") < 3 x 10'". 

PACS numbers: 12.15.Ji, 13.25.Hw, 14.40.Nd 



We report results of searches for B decays to the 
charmless two-body final states Q wtt+j i?+ 

ujK+, B+ r/'7r+, B+ t]'K+, B° ujK°, B° W7r°, 
and rj'K'^ . These processes are manifestations of 

penguin or suppressed tree amplitudes proportional to 
small couplings in hadronic flavor mixing (CKM matrix 
0|). Because of the absence of CKM favored b ^ c am- 
plitudes these decays are particularly sensitive to poten- 
tially new contributions from interference effects and vir- 
tual particles in loops. Previous measurements yielded 
an unexpectedly large rate for B rj'K, motivating a 
number of new theoretical ideas. The precise measure- 
ment of these and additional rare B decay modes will 
enable a better understanding of the underlying decay 
mechanism, including the possible contribution of physics 
beyond the standard model. This in turn will contribute 
to the measurement of fundamental parameters, includ- 
ing the CP-violating CKM phases. 

The data were collected with the BABAR detector Q 
at the PEP-II asymmetric e+e' collider 5] located at 
the Stanford Linear Accelerator Center. The results pre- 
sented in this paper are based on data taken in the 1999- 
2000 run. An integrated luminosity of 20.7 fb'^, corre- 
sponding to 22.7 million BB pairs, was recorded at the 
T{AS) resonance ( "on-resonance" , 10.58 GeV), with an 
additional 2.6 fb'^ about 40 MeV below this energy ( "off- 
resonance") for the study of continuum backgrounds. 

The asymmetric beam configuration in the laboratory 
frame provides a boost to the T{AS) increasing the mo- 
mentum range of the i3-meson decay products up to 
4.3 GeV/c. Charged particles are detected and their 
momenta measured by a combination of a silicon vertex 
tracker (SVT), consisting of five layers of double-sided 
detectors, and a 40-layer central drift chamber (DCH), 
both operating in the 1.5 T magnetic field of a solenoid. 
Photons and electrons are detected by a CsI(Tl) electro- 
magnetic calorimeter (EMC), which provides excellent 
angular and energy resolution with high efficiency for en- 
ergies above 20 MeV 

Charged particle identification (PID) is provided by 
the average energy loss (dE/dx) in the tracking de- 
vices and by a unique, internally reflecting ring imaging 
Cherenkov detector (DIRC) covering the central region. 
A Cherenkov angle K-tt separation of better than 4 stan- 
dard deviations {a) is achieved for tracks below 3 GeV/c 
momentum, decreasing to 2.5a at the highest momenta 
in the final states considered here . Electrons are iden- 
tified with the use of the EMC. 

We reconstruct a B meson candidate by combining an 



uj or T]' candidate with a charged track, tt*^ — ^ 77, or 
A'g TT+TT'. The resonance decays R we reconstruct 
are u; 7r+7r-7r°, 77' 7]tt+tt- (77^,^^), or 'n'^p°-f {r]'p^), 
with rj ^ and p'^ — > tt^tt'. These modes are kine- 
matically distinct from the dominant B decays to heav- 
ier charmed states. Backgrounds come primarily from 
combinatorics among continuum events in which a light 
quark pair is produced instead of an T(4S'). 

Monte Carlo (MC) simulations of the target decay 
modes and of continuum background are used to establish 
the event selection criteria. The selection is designed to 
achieve high efficiency and retain sidebands sufficient to 
characterize the background for subsequent fitting. Pho- 
tons must satisfy > 50 (100) MeV for tt° (r?) can- 
didates. For rj'^p^^ candidates from i?"*" -q' and 
5+ rj'iT^ the requirement is E^ > 200 MeV, while 
for B° ri'K° it is looser {E^ > 100 MeV) because of 
smaller combinatoric background. 

We select uj, rj', rj, and p candidates with the follow- 
ing requirements on the invariant masses in MeV/c^ of 
their final states: 735 < TO(7r+7r-7r°) < 830, 930 < 
TO(r77r+7r-) < 990, 900 < m{p-f) < 1000, 490 < 771(77) < 
600, and 500 < m{T:+7:-) < 995. For vr" and K'^ 
candidates we require 120 < m{j'-f) < 150 and 488 < 
m(7r+7r-) < 508. 

Tracks in lo, rj', or p candidates must have DIRC, 
dE/dx, and EMC responses consistent with pions. For 
charged B decays, the B primary track must have an 
associated DIRC Cherenkov angle within 3.5(7 of the ex- 
pected value for a kaon or pion. For modes with 
the three-dimensional fiight distance from the produc- 
tion point must exceed 2 mm, and the angle between the 
fiight and momentum vectors projected perpendicular to 
the beam must be less than 40 mrad. 

A B meson candidate is characterized by two kinematic 
observables. The minimally correlated pair we use are 
the energy constrained mass 77t,ec and energy difference 
AE. In the T{AS) frame the B meson energy E* equals 
the beam energy E^^^^^. A kinematic fit of the measured 
candidate four momentum in this frame with the con- 
straint E* = -E^ea,^ yields mEc, while AE = E* - E^^^^ 
measures the consistency of this constraint. We require 
\AE\ < 0.2 GeV, and niEC > 5.2 GeV/c^. The resolu- 
tions on these quantities are mode dependent but average 
about 30 MeV and 2.8 MeV/c^, respectively. 

To discriminate against tau-pair and two-photon back- 
ground we require the event to contain at least five 
charged tracks. To reject continuum background we 
make use of the angle 9t between the thrust axes of the B 
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candidate and the rest of the tracks and neutral clusters 
in the event, calculated in the center-of-mass frame. The 
distribution of cos 9t is sharply peaked near ±1 for com- 
binations drawn from jetlike qq pairs, and nearly uniform 
for the isotropic B meson decays. 

The yields are obtained from extended unbinncd max- 
imum likelihood (ML) fits, with two variants specified in 
the following paragraphs. The first (MLl), which pro- 
vides our results for all modes except 5° ujtt'-', uses 
several uncorrelated variables for the kinematics of the 
B decay chain and a Fisher discriminant for the produc- 
tion and energy flow. The second (ML2) is applied to all 
channels with an to meson; it uses AE and the output of 
a neural network built from the remaining inputs. Com- 
parisons for the ujTT^ , ujK^ , and ujK^ modes show that 
the central values and errors for the yields obtained by 
the two approaches are in very good agreement. Simple 
cut-based analyses are performed as checks for each final 
state. Agreement of central values is good in all cases, 
although, as expected, errors are larger than for the ML 
analyses, particularly for modes having high background. 

The MLl fit method is applied to events satisfying 
\cos9t\ < 0.9. The input observables are Ai?, toec, 
the invariant mass mn of the intermediate resonance, the 
Fisher discriminant T, and, where relevant, the rj mass 
TO,;, the measured DIRC Cherenkov angle for the B pri- 
mary track, and the cosine Ti of the helicity angle, the 
angle in the lo rest frame between the normal to the lo 
decay plane and the B flight direction. The Fisher dis- 
criminant combines eleven variables: the angles with 
respect to the beam axis in the T{AS) frame of the B 
momentum and B two-body decay axis, and a nine bin 
representation of the energy flow about the B decay axis. 

For the ML2 method we relax the preliminary require- 
ments to 100 < m(77) < 160 MeV/c^ and |A£:| < 0.3 
GeV. The neural network is constructed with the B mo- 
mentum p*, a for resonance masses, 7i, and variables 
representing energy flow and angular distributions, in- 
cluding 6t- 

We use MC to estimate backgrounds from other B de- 
cays, including final states with and without charm. For 
most of our modes we find contributions that are negli- 
gible. For the rj'—^p'^^ modes we account for small cross 
feed contributions in the systematic error estimate. 

The likelihood function for N events is 



TV! 



Here Uj is the population size for species j (e.g., signal, 
background) and Vj{:x.i) the corresponding probability 
distribution function (PDF), evaluated with the observ- 
ables Xi of the ith event. 

For the fits of charged B decays Ci becomes 

nc [fKcVKci^i] + (1 - /kc)Kc(xj)] , 



where nT^{nK) is the number of _B+ Rtt^{B^ RK^) 
signal events, nc is the number of continuum background 
events, and fKC is the fraction of continuum background 
events for which the B primary track is identified as a 
kaon. These quantities are the free parameters of the 
ML fit. The probabilities for the components are V-j^s 
{Vks) for B+ R'K+{B+ RK+) signal and V^c 
(Vkc) for background where the primary track is a pion 
(kaon) . Since we measure the correlations among the ob- 
servables in the data to be small, we take each Vj to 
be a product of the PDFs for the separate observables. 
The analyses involving a K'^ are treated identically ex- 
cept that there is only one component of signal and of 
continuum background. 

A second B candidate satisfying the preliminary cuts 
occurs in about 10-20% of the events. In this case the 
"best" combination is selected according to a quantity 
computed with toec, "t-Aj nT-rj (foi' v' ~* rjir^Tr^ modes), 
and the Fisher discriminant. 
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FIG. 1: Invariant mass distributions for inclusive data sam- 
ples of candidates with T(4S') frame momentum greater than 
2.3 GeV/c for (a) 77', with 520 < m(77) < 575 MeV/c^ and 
(b) Lu candidates, with 120 < 711(77) < 150 MeV/c^. From 
the overlaid fit curves the Gaussian peak widths are 4 and 
10 MeV/c^, respectively. 

We determine the PDFs for the likelihood fits from 
simulation for the signal component, and from off- 
resonance and sideband data for the continuum back- 
ground. Peaking distributions (signal masses, AE, J-) 
are parameterized as Gaussians, with or without a second 
Gaussian or asymmetric width as required to describe 
the distributions. Slowly varying distributions (combi- 
natoric background under mass or energy peaks, 7i, or 
JF) have first or second order polynomial shapes. The 
combinatoric background in toec is described by a phase 
space motivated empirical function j^. Control samples 
of B decays to charmed final states of similar topology 
are used to verify the simulated resolutions in AE and 
Toec- Inclusive resonance production samples such as 



those shown in Fig. ^ are used similarly for the relevant 
B daughter mass spectra. 



TABLE I: Signal event yield with statistical uncertainty, de- 
tection efficiency e, daughter branching fractions that were 
forced to 100% in our signal mode simulations, significance S 
(defined in the text), and branching fraction result for each 
decay chain or mode, with the final (combined) result given 
in bold type. We show 90% confidence upper limits in paren- 
theses where appropriate. 



Mode 



Yield 



e UBi S 

% % a 



10' 



rj'K+ 



49.5 
87.6 



8.1 
7.3 
+ 13.4 
12.5 



20 17.4 15 
18 29.5 11 
17 



63^ 



-10 

80f^ 



-11 

70 ± 8 ± 5 



20.8+_l:t 



16 
16 



6.0 4.7 
10.1 4.2 
5.9 



j + 15 



28: 



42 



-1-13 
-11 



19 



±4 



5.7 



-f3.8 



-2.8 



20 
19 



17.4 3.2 

29.5 0.1 
2.8 



7.1 



-1-4.8 
-3.5 



5.4+11 ± 0.8 (< 12) 



6.4 



-1-5.6 
4.4 
1-1-4.6 



27.61?:? 
-0.9 



-fS.O 
3.2 



22 88.8 1.3 

18 30.5 3.2 

21 88.8 4.9 

18 88.8 



1.4+l;l ± 0.3 (< 4) 
6.4^2 8 ± 0.8 (< 13) 
g ± 0.7 
-0.3 ± 1.1 ± 0.3 (< 3) 



We compute the branching fractions from the fitted 
signal event yields, reconstruction efficiency, daughter 
branching fractions, and the number of produced B 
mesons, assuming equal production rates of charged and 
neutral pairs. To determine the reconstruction efficiency, 
including any yield bias of the likelihood fit, we apply the 
method to simulated samples with the signal and contin- 
uum background populations expected in the data. Ta- 
ble n shows for each decay chain the branching fraction 
we measure, together with the quantities entering into 
its computation. The statistical error on the number of 
events is taken as the shift from the central value that 
changes the quantity = — 2 log (£/£inax) by one unit. 
We also give the significance S, computed as the square 
root of the difference between the value of for zero 
signal and the value at its minimum. The used for 
significance includes a term that accounts for the addi- 
tive systematic error discussed below. Where the signifi- 
cance is less than four standard deviations, we quote also 
(Bayesian) 90% C.L. upper limits, defined by the solution 
B to the condition C{b)db/ C{b)db 0.9. 

In Fig. 121 we show projections of toec and AE for 
the modes with significant yields. The projections are 
made by selecting events with signal likelihood (com- 
puted without the variable plotted) exceeding a mode- 
dependent threshold that optimizes the expected sensi- 
tivity. 
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FIG. 2: B candidate rriEC and AE for B+ ??'-R'+(a, b), 
i3" Tj' K'''{c, d), and ti;7r"'"(e, f). Histograms represent 

data, with the 77' ^ rirnr subset shaded, solid curves repre- 
sent the full fit functions, and dashed curves the background 
functions. 



We have evaluated systematic errors, which are dom- 
inated in most cases by the PDF uncertainties (3-18%, 
depending on the decay mode). To determine these we 
vary parameters of the PDFs within their uncertainties 
and observe the impact on the fit yield. We include them 
in upper limits by convolution with the likelihood func- 
tion. This is the only additive systematic error; all others 
are multiplicative. The estimate of any systematic bias 
from the fitter itself (1-4%) comes from fits of simulated 
samples with varying background populations. 

The uncertainty in our knowledge of the efficiency is 
found from auxiliary studies to be 1% per track, 1.25% 
per photon, and 5% per K*^ for the candidate B and the 
unreconstructed B, which must contribute tracks to fulfill 
the event multiplicity requirement. We add these errors 
linearly for the required tracks in the event, and similarly 
for the photons and neutral kaons. Our estimate of the 
B production systematic error is 1.6%. Published world 
averages |lOl| provide the B daughter branching fraction 
uncertainties. 

Systematic errors associated with the event selection 
are minimal given the generally loose requirements. We 
account explicitly for cos^t (1%), for which we observe a 
nearly uniform distribution in the signal simulation. We 
also include errors of 4% from those PID requirements 
that are imposed via cuts rather than the fit. 

We have observed signals of at least 4(7 in five of the 
decay chains studied here, as reported in TableQ] Where 
we have multiple chains for a given mode we combine 
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the results by adding the distributions that represent 
them and their uncorrelated statistical and systematic 
errors. 

The final results are generally in agreement with those 
previously reported |3llll|. with somewhat smaller errors. 
In particular, we confirm the expected B{B^ a;7r+}> 
B{B~^ ujK^), and the rather larger than predicted 
rate for B — > rj'K obtained by the CLEO Collaboration 
. Coiyectured sources of ij' enhancement include fiavor 
singlet charm enhanced and constructively in- 
terfering internal penguin diagrams [l^ Il5l | . Our results 
in combination with expected measurements of related 
modes involving rj and K* should help to clarify this sit- 
uation. 
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